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Wet sand flows better than dry sand 
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We investigated the yield stress and the apparent viscosity of sand with and without small amounts 
of liquid. By pushing the sand through a tube with an enforced Poiseuille like profile we minimize 
the effect of avalanches and shear localization. We find that the system starts to flow when a critical 
shear of the order of one particle diameter is exceeded. In contrast to common believe, we observe 
that the resistance against the flow of wet sand is much smaller than that of dry sand. For the 
dissipative flow we propose a non-equilibrium state equation for granular fluids. 

PACS numbers: 05.70.Ln, 83.60.La, 45.70.Mg, 83.80.Fg. 



Capillary forces between single grains are the cause of 
the stiffness of sculptured wet sand, as opposed to dry 
sand which does not support its own weight There- 
fore, it is commonly believed that wet sand should be 
more resistant, i.e. more "viscous", than dry sand. Not 
many studies, though, on the flow behavior of wet sand 
have been undertaken yet. The few studies at hand have 
been performed mostly either in a rotating drum Q or in 
standard rheometric set-ups Our study was inspired 
by an apparatus invented by the authors of references 
PL H @ > but their work is more focused on the regime be- 
low the yield point. Besides its practical applications in 
e.g. civil engineering or industrial processes, the Theolog- 
ical or dynamical properties of dense packed granular ma- 
terials are investigated to help develop a non-equilibrium 
thermodynamic formalism. Whilst the description of a 
rapid flowing granular gas is well established [7[, the op- 
posite limit of a dense packing remains an open ques- 
tion. Proposals for state equations for granular matter 
include, for example, the Enskog equation either for de- 
scribing rapid dilute granular flows Q or for the high 
density limit [9( . Yet, it is not possible to describe all 
of the experimental findings, e.g. the fact th at g ranular 
matter does not follow a Boltzmann statistic [l(J HH • 

Granular matter is an example of a so-called yield 
stress fluid. Above a critical yield stress, the mechani- 
cal response changes from a solid-like or jammed state 
to a flowing one. An illustrative study on the rheolog- 
ical prop erties of wet soils was given, for example, in 
ref. [jjj]- These authors measured viscoelastic proper- 
ties under steady and oscillatory stress. Soils are a typi- 
cal example of wet granular material, and they exhibit a 
mixed behavior of elasticity, viscosity and plasticity that 
is mainly determined by the quantity of water, the yield 
stress, shear modulus, and viscosity decrease with the wa- 
ter content. The water content is changed by adding liq- 
uid to a given sample, thus reducing the packing density 
of the solid grains. The system starts to flow when the 
external stress exceeds the inter-aggregate contact [l3[ 
forces. The mechanical properties at lower water content 



'Electronic address: c.wagner@mx.uni-saarland.de 



are determined by the liquid bridges between clay clus- 
ters and quasi-crystals, and those at higher water content 
are determined by the flow of the liquid through the soil 
pores (p3 |. 

A principle difficulty in the determination of the yield 
stress point arises from the fact that a yield stress fluids 
exhibits thixotropy, also referred to as aging or rejuve- 
nation [|| . Yield stress and thixotropic behavior depend 
on the degree of compaction and the topology. Both 
can be modulated by weak attractive forces like the elec- 
trostatic charge and humidity [ijj]- Macroscopically, the 
jammed nature of dry and wet granular materials de- 
termines the viscoelastic properties, e.g. the route to 
the jammed transition can be explained by a modified 
Vogel-Fulcher-Tammann behavior [l5|, [l6| , and it can be 
understood by the relation between the compaction and 
the grain mobility (T7L [lH . In dense granular materials 
the resistance to flow and the dissipation of energy oc- 
curs when the bulk granulate has to pass from one con- 
figuration to another. This rearrangement process leads 
to a new, jammed configuration. The rearrangement dy- 
namics driven by shearing causes the rejuvenation and 
aging of the fluid-structure and it is meaningful to de- 
termine the yield stress point when the structural in- 
tegrity is maintained and shear localization, avalanches, 
and arching effects are prevented. We will show that by 
doing so, it is possible to reach a steady state with a con- 
stant relation between the pumping and dissipation of 
energy if the granular material is sheared in a oscillatory 
way. Our sand was composed of glass spheres of diameter 
D = 140 — 150/im with and without additional dcionized 
water. The content of water in the paste W is defined as 
the quotient between the liquid volume and the volume 
occupied by the grains. As a first attempt to measure the 
yield stress point of our samples, we followed the meth- 
ods used to study the mechanics of soil. Measurements 
were done with a rotational viscometer (MARS, Thermo 
Haake. Karlsruhe, Germany) with a parallel plate geome- 
try 1 1 21 ] , and with sand paper glued on the on the surfaces 
of the plates. Further improvement was given by using 
a vane sensor or a "bars" sensor, which is shown in the 
inset of fig. [TJ However, for our study, all geometries 
gave comparable results, while the latter was the most 
reproducible. A typical protocol to study yield stress 
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FIG. 1: Comparative rheological test under constant defor- 
mation (CD modus) for micron-sized dry sand (filled black 
squares) and the wet paste with a water content of W = 0.05 
(open triangles). In the inset, the sensor for the "bars" ge- 
ometry is shown. 



fluids - for example coatings and paints - involves im- 
posing a constant shear rate and measuring the strain. 
The so-called yield stress point is the highest stress at 
which no flow occurs at a given shear rate, e.g., at the 
maximum in the flow curve of fig. [TJ Beyond the max- 
imum, the deformation is large enough to compromise 
the structural integrity of the granular fluid. Figure 1 
shows a representative result for the case studies of dry 
and wet sand. The maximum corresponds to the stress 
at which the material begins to plastically deform; be- 
fore that, the material deforms elastically. The test must 
be carefully tuned, because the shear rate must be small 
enough to catch the elastic and plastic regimes but, if it is 
too small, the peak is difficult to identify. In general, the 
yield stress point depends on the time scales and shear 
rates applied in the experiment - however, it was recently 
found that the critical yield strain of many different yield 
stress fluids, like ketchup or pastes, does not depend on 
the external parameters; instead, it depends only on the 
type of fluid [H. 

The yield stress point for the dry sand is t = 0.8fcPa, 
7 C = 0.1 = D/10 (r is the stress and 7 the strain) while, 
for the wet sand, t = 4.23/cPa, 7 C = 1.4 = 3D/2, for a 
given shear rate of 7 = 0.0038 s _1 and 7 = 0.003 s _1 
respectively. Apparently, the liquid bridges between the 
grains reinforce the structural integrity of the wet sand. 
Similarly, one needs less torque to turn the sensor at 
a given rotational speed in the dry sand than the wet 
sand. The asymptotic values of the stress are 0.2 and 
1/cPa respectively; however, due to the breakdown of the 
microstructure and flow localization, one cannot straight- 
forwardly deduce the yield stress (or even a shear viscos- 
ity) later on. To override this problem, we carried out ex- 
periments with the shear cell proposed in ref. 0, HI, @ (see 
the inset of fig. [2]). The use of this device made it possi- 
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FIG. 2: Family of differential pressure curves for increasing 
shear amplitudes (see text for further explanation); the inset 
shows a schema of the measurement cell. 



ble to shear the pastes all over the sample volume. The 
granular material is put into a acrylic glass cylindrical 
cell at atmospheric pressure. The sides of the cell are 
sealed with thin latex membranes of 300/im, which are 
flat at the beginning of the experiment. Chambers filled 
with water are attached to both sides of the cell, ad- 
jacent to the membranes, where it is possible to inject 
and extract water into and out of the chamber through 
tubes via a syringe. The syringe-pistons are connected 
to spindles that could be moved with a stepper motor. 
The pistons move at the same speed as each other but in 
the opposite direction to prevent dilatation. The acrylic 
glass tube we used for this experiment was 24 mm long. 
It was filled with dry sand or with a wet sand paste with 
W = 0.03 so that, in both cases, the global packing frac- 
tion pf = 0.63. We should mention that this number 
must be considered cautiously, since we ignore the local 
compaction of the granulate, e.g., the topology of poly- 
tetrahedral grain structures [2fj |. 

The deformation of the membranes is approximately 
parabolic and imposes a Poiseuille-like flow profile within 
the sample; thus, avalanching and arching can be kept 
to a minimum. Piezoresistive sensors sense the pressure 
(with respect to atmospheric pressure) in both cham- 
bers. A fixed quantity of water AT^ is injected or ex- 
tracted in/from the adjacent chambers, and the differ- 
ential pressure p\ — pi is measured (see the inset of fig. 
[2]). The displacement of the membrane is calculated as 
Ax = AAV/ (tt-Dc), where D c is the diameter of the 
acrylic glass cell (24mm). The paste is sheared at a 
very low shear rate with a shear amplitude a s , between 
—a s and a s . The protocol was to move the membrane 
Ax = 6/iTO in one second, wait 13 seconds, and measure 
Pi and pi until the selected a s is reached. A curve of the 
pressure difference p\ — pi as a function of the displace- 
ment - which represents the strain - is obtained for each 



3 



□_ 




* -5 



51 

5 
-10 



-1000 -500 500 1000 

(c) 




-50 







50 



40 


-401 




-1000 -500 500 1000 

x (jam) 



FIG. 3: Differential pressure versus the shear displacement. 
(a),(b): A sample of 17.7 gr. of dry sand (26% room 
humidity,26C°) is compared with (c),(d): a paste of the same 
sand content mixed with water (W = 0.03). 



shearing amplitude a s . A family of differential pressure 
characteristics for increasing shear amplitudes are shown 
in fig. [5] The slope of the two branches corresponds to 
the restoring force exerted by the rubber membranes, but 
a stress r s t is stored in the system which is related to the 
opening of the loop with 2T s t(a s ) = A(pi — p2)o- In fig. 
[31 the differential pressure curves for dry and wet sand 
are shown for the two limiting cases of small and large 
shear amplitudes. By shearing with a s < 2D, we find 
that the dry sand does not resist any stress (fig. El (a)), 
while the wet sand behaves like a yield stress fluid due 
to the liquid bridge network (LBN) (fig. [3J (c) and the 
inset of fig. [3]), with a A(pi — p2)o ~ IfcPa. Following 
ref.[3, 0, 0, the capillary pressure is an energy density 
p = dE/dV < pfNj/D where N = 6 is the LBN coor- 
dination number, 7 is the surface tension of water, and 
D = 145um; thus, p < 1.9/cPa and, in agreement with 
we find that the restoring forces originate 
from the capillary forces between the grains. The onset 
of dissipative flow occurs at a s > 2D then the opening 
of the loop starts to increase. It grows more strongly 
for the dry than for the wet sand, with maximum values 
of 116/cPa, and 22/cPa for the largest shear amplitude, 
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FIG. 4: Stress - strain behavior for increasing shear amplitude 
for dry (triangles) and wet sand W = 0.03 (squares) and the 
corresponding fittings for the dissipative flow range (equation 
Q). The inset shows, for the wet sand, the range where the 
stored stress is the yield stress. 



respectively. This means that one needs less energy to 
push the wet sand through the tube. It is worth men- 
tioning that, for a given shear amplitude, the hysteresis 
loop is stable during a many cycles. We carried out a 
few tests which were one week in duration and found 
that the opening of the loop for different shear ampli- 
tudes remains constant (we made sure that all our data 
were taken in such a steady state situation). We also 
tested shear velocities between 0.6 and 300/im/s without 
observing differences in the hysteresis loop. 

Figure 3] shows the complete stress-strain behavior for 
the dry and the wet sand. The inset shows the differential 
pressure for wet sand for a s < 2D, from which we deduce 
a mean value for the yield stress r = (A(pi — ^2)0} /2 = 
(0.6 ± 0.5) fcPa. For the dry sand, we find (in this range) 
(A(pi — p2)o) = (0 ± 1) fcPa. In the dissipative flow 
range a s > 2D, it is not possible to shear the sample 
without changing the local compaction by rearrangement 
of the grains, and the response of the system is a change 
of configurations between jammed states. The experi- 
ment shows that the addition of a small amount of wa- 
ter changes the flow resistance dramatically, possibly due 
to the existence of lubrication layers and because it fa- 
vors cluster formation that minimizes the effective free 
volume [20]. Near the onset of yield, the stress-strain 
curve for increasing a s (fig. f?]) can be modeled with a 
stretched exponential: 
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A fit with with equation |T]) yields an onset value D* = 
150/im for both cases which is close to the mean diameter 
of the grains; for the constant Cq and the exponent n, we 
find cq = 55, n = 1 and Cq = lOfcPa, n — 3/2, for the 
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FIG. 5: Stress - strain curve for wet sand obtained by in- 
creasing (triangles up) and decreasing the shear amplitude 
(triangles down). 



dry and wet sand respectively. The work that is applied 
to the system is given by W = F.a s , and its response 
is related to the energy E that is dissipated during the 
rearrangement mechanism; thus, it is possible to rewrite 
cq. H|) as 



W_ 

Wo 



exp 



E 



(2) 



where Wo = 2F.D* and E$ is the activation energy 
for the rearrangement. The latter should be related to a 
configurational entropy; we find Eq = 88 and 16nJ per 
grain for the dry and wet granulate respectively. 

For n — 1, eq. {I) it is a Doolittle-like equation [2l|. It 
was shown that the viscosity of a granular fluid depends 
strongly on the particle concentration [22I ] and, in terms 
of a free volume theory, E oc 1/vf. The activation en- 
ergy is related to the global packing factor Eq oc pf (the 
occupied volume divided by the cell volume). Following 
emerging concepts on non-equilibrium thermodynamics 
[HI, we propose that the equation Q could be a state 



equation for granular fluids in dissipative flow under glob- 
ally constant packing fraction conditions. 

It is difficult to estimate the value of Eq for dry sand 
since the energy loss is sensitive to, e.g., microscopic 
roughness and topology. But, for the wet sand, it can 
be related to the liquid bridges. By considering the co- 
ordination number N — 6, the dissipation per bridge is 
E /N = (2.6 ± 0.5)nJ, which is in good agreement with 
the estimation for the energy loss during bridge forma- 
tion and rupture AE cap < ■K"fD 2 /2 = 2AnJ 6]. Figure 
[5] shows a stress-strain curve for the first scan of increas- 
ing and decreasing the shear amplitude applied to wet 
sand. Beyond a s > 2D, the structure of the granu- 
lar fluid changes during the successive rearrangements. 
Possible mechanisms for this include the formation of 
water channels by coalescence of the liquid bridges and 
larger polytetreahedral solid structures. The decreasing 
branch of Fig. 5 might thus correspond to a fluid that ex- 
hibits phase separation and shear bands. For the range 
a s < 2D, the opening of the loop (the yield stress) tends 
to zero; again, this is an indication for a phase separa- 
tion. Similar measurements for the dry sand revealed no 
difference between increasing and decreasing shear am- 
plitude, but the rearrangement process there occurs not 
between evolving clusters, but between the single grains. 

In conclusion, we have found that it is much easier to 
push wet sand than dry granular matter in a Poiscuillc- 
like profile through a tube. Even if the capillary forces in- 
crease the yield stresses, it looks like the water promotes 
the cluster formation and reduces the inter-grain friction. 
The leading dissipation mechanism results from the rup- 
ture and formation of liquid bridges, and we are able to 
explain our data quantitatively within the framework of 
an excluded volume theory. Finally, we find indications 
that the yield of the system is related to the microscopic 
size of the grains. Further studies with different grain 
sizes are ongoing. 
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